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Abstract 



We calculate the super symmetric 0(a s ) QCD corrections to stop and sbottom decays 
into vector bosons within the Minimal Supersymmetric Standard Model. We give analytic 
formulae and perform a numerical analysis of these decays. We find that SUSY-QCD 
corrections to the decay widths are typically —5% to —10% depending on the squark 
masses, squark mixing angles, and the gluino mass. 



1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) [Q every quark has two scalar part- 
ners, the squarks q L and q R . In general q L and q R mix to form mass eigenstates q ± and q 2 (with 
< rriq 2 ), the size of the mixing being proportional to the mass of the corresponding quark 
q 0. Therefore, the scalar partners of the top quark (stops) are expected to be strongly mixed 
so that one mass eigenstate ti can be rather light and the other one t 2 heavy. The sbottoms 
bi and b R may also considerably mix for large tan/3 = v 2 /v% (where v% and v 2 are the vacuum 
expectation values of the two Higgs doublets). 

Squark pair production in e + e~ annihilation including mixing has been studied at tree-level in 
Ref. ||, then including conventional QCD corrections in j|, [5], and supersymmetric (SUSY) 
QCD corrections in Ref. |7j. The cross section of squark production at hadron colliders (LHC 
and Tevatron) in next-to-leading order of SUSY-QCD was given in Ref. ||. 

While, quite naturally, most studies concentrated on t\ (b\) production and decays, those of 



the heavier t 2 (b 2 ) have been discussed much less || [H], [11]. These particles could be produced 
at the LHC or an e + e~ Linear Collider. The decay patterns of the heavier mass eigenstates t 2 
(b 2 ) can be very complex due to the many possible open decay channels. There are the decays 
(i,j = l,2;k = l..A): 

ti^tx° k ,bxj, k -^bxt txj, (!) 
and in case the mass splitting is large enough: 



t i ^b j W+, bi^tjW, 

t 2 -> i x z°, ~b 2 -> hi z°, 

and 

t 2 -> h (h°, H°, A ), b 2 ^ h (h°, H°, A ). 



(2) 
(3) 



Stops and sbottoms can also decay strongly through ti — > tg and b\ bg. If these decays are 
kinematically possible, they are important. The SUSY-QCD corrections to squark decays into 
charginos and neutralinos have been calculated in |Q. The SUSY-QCD corrections to squark 



decays into Higgs bosons have been treated in [13|, and those into gluino in [f[4 |. The decays 
to photon and gluon, q 2 — > qij, qig, which are absent at tree-level, are not induced by 0(a s ) 
SUSY-QCD corrections either. This is due to SU(3) C x U(l) em gauge invariance. 

In this paper we calculate the 0(a s ) SUSY-QCD corrections to the squark decays into W ± and 
Z° bosons of eq. @. The stop decays into vector bosons can be dominant in the parameter 
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region where (i) rrit (At — fi cot/3), appearing in the mass matrix, is large enough to give the 
necessary mass splitting, (ii) M and are relatively large to suppress the decay modes to 
chargino, neutralino, and gluino, and (iii) the decays into Higgs particles are not so important 
(for instance, itla is large). For sbottoms to decay into vector bosons instead of (i) it is necessary 
that mb (A b — fM tan/3) is large enough to lead to a large mass splitting of b\ and b 2 - 

In calculating the SUSY-QCD corrections we work in the on-shell renormalization scheme and 
use dimensional reduction []15|| , which preserves supersymmetry at least at two-loop level. The 
gluon loop corrections are calculated in the Feynman gauge. The couplings of the squarks to 
vector bosons depend on the squark mixing angles for which an appropriate renormalization 
procedure is necessary. Here we use the renormalization prescription as introduced in |7| , where 
it was applied to the case of e + e~ — > q^,. 



2 Tree— level formulae and notation 



The current eigenstates q L and q R are related to their mass eigenstates q 1 and q 2 by 

(*}\=n*(*A K 9 =( COS6 " Sin M (O<0 ff <7r). (4) 

In the (q l3 q 2 ) basis the squark interactions with Z° and W bosons are given by: 

C = ^-Z^KiKj^ + ^K^C^qj&q, 

cos Ow 

n\ x n) x w- b) # t, - ^= n\{R} n w+ 1] ? K (5) 

where C^r = J| L R — e q sin 2 W with J| the third component of the weak isospin and e g the 
electric charge. Here it is assumed that the (3,3) element of the super-CKM matrix ^33 = 1. 

At tree-level the amplitude of a squark decay into a W or Z° boson has the general form 

M (q? - tfjV) = -igc l3V (h + hre;(k 3 ), (6) 

with ki, k 2 , and k 3 the four-momenta of qf, q?, and the vector boson V (V = W ± ,Z°), 
respectively (Fig. la), a and (3 are flavor indices. In the following we define = m^, 

~a ~(3 

rrij = m~0 , TLq. = 1Z ik , TZjk = 7Z jk for simplicity. Moreover, we shall use primes to explicitly 
distinguish between different flavors. 
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With this notation the qf-q^-V couplings Cyy are: 



i ( IIl c °s 2 9q ~ e <? sin2 °w -§ IL sin 2 % 
cos e ^ ^ - 1 Jf L S in 2fl 9 ~ J| L sin 2 ^ - e g sin 2 ^ w 



ij 



_ 9 2 (c ijV ) 2 K 3 (m 2 ,m 2 ,rr#) 
167rr?Ty to? 



(7) 



- J_-7? 7?' - if cos0 ? cos0y - cos 0$ sin 0y \ 
( i( „ ^'UiKji - j5 [ - s in^-cos^ sin ^ sin ^ ) ' l ' 



The tree-level decay width can thus be written as 



with k(x, y, z) = (x 2 + y 2 + z 2 — 2xy — 2x2 — 2yz) l l 2 . 

Note that the decay q L m\ — >■ <?_r(l) ^° does n °t occur at tree-level, and q { — > is not 

allowed at tree-level if one of the squarks is a pure "right" state. 



3 SUSY-QCD corrections 

The 0(a s ) loop corrected decay amplitude is obtained by the shift: 

CijV -> Cyy + 5<4j(r + 5cjjy + 5 %°y> ( 10 ) 

in eq. (0). The superscript v denotes the vertex correction (Figs, lb-f), w the wave function 
correction (Figs, lg-i), and c the counterterm to the couplings. Thus the 0(a s ) corrected 
decay width T is given by 

r = r° + <yrw + sr^ + 5r (c) + 5r rea/ (11) 

where 

, v q 2 K 3 (rrif, m 2 , rnh) r r„\ > 

87T TOy - 1 

5r reaZ is the correction due to real gluon emission (Figs, lj-^) which is included in order to cancel 
the infrared divergence. The total correction can be decomposed into a gluon (exchange and 
emission) contribution 8T^ 9 \ and a gluino-exchange contribution 5T^': T = T° + 6T^ 9 ' + 5T^\ 
The contribution from squark loops (Figs. 1 f,i) cancels out in our renormalization scheme. 
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3.1 Vertex corrections 



The vertex corrections stem from the five diagrams shown in Figs. 1 b-f. The gluon-exchange 
graphs of Figs, lb-d yield 



0C ijV 



3^ 



B (m 2 , A 2 , m 2 ) + B (m 2 , A 2 , m 2 ) 



- 2 (m 2 + m 2 - m 2 ) (C + d + C 2 ) 



(13) 



i?o and 6*0,1,2 are the standard two- and three-point functions for which we follow the 
conventions of |T7|]. Here we use the abbreviation C m = C m (m 2 ,my,'m 2 ] A 2 ,m 2 ,m 2 ). A gluon 
mass A is introduced to regularize the infrared divergence. 

The gluino-exchange contribution, Fig. le, gives 



Sc. 



21Z 



( fh ^ {l% L [ 2 m 2 C + m 2 2 d + m\ C 2 + (m 2 - m 2 ) (d + C 2 ) + B ] sin 2^- 

(14) 



+ 2m~ g m q {Il L - 2e q sin 2 W ) (C + d + C 2 ) cos 2^}, 



with C m = C m (m? ml,m\-m\,m\ 1 m 2 ) and 5 = B (m^,m 2 m 2 ), for the decay g 2 -> g x Z c 



and 



Sc. 



ijW 



O'.s 



v/2 



3 7T 



{m 5 (C + Ci + C 2 ) (m, ft i2 72^ + nv 7^ 2 ) 

- [ m 2 d + m 2 d + m 2 (2C + d + d ) + B ] K a % x 

- m q raj (d + d) ^ i2 7^- 2 } 



(15) 



with C m = C m (m 2 , 777^, m 2 ; m|, m 2 , m 2 ,) and £? = B (m^, m 2 , m 2 ,), for the decay ^ — > H /± . 
The squark loop of Fig. If does not contribute because it is proportional to the four-momentum 
of the vector boson. The total vertex correction is thus given by: 



0C ijV 



0C ijV + 0C ijV ■ 



(16) 



3.2 Wave— function correction 



The wave-function correction is given by 



Sc 



ijV 



SZ li (q^)+SZ : 



c ijV + SZik(q?)c kj v + SZji(q?)cav, (i^kj^l). (17) 
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Z n m{<In) are the squark wave-function renormalization constants. They stem from gluon, gluino, 
and squark loops (Figs. 1 g-i) and are given by: 

Re Is^'^fm 2 )\ 

8Z™m = -Ite 6Z%p(q n ) = - L_S \ (18) 

flux ''^n , 

with S„ n (m 2 ) = dT lnn (p 2 )/dp 2 \ p 2 =m 2. (The gluon loop due to the qqgg interaction gives no 
contribution because it is proportional to the gluon mass A — > 0.) The squark self-energy 
contribution due to gluon exchange (Fig. lg) is 



2as 

3 7T 



5 (m| n , 0, m 2 J + 2m| n S (m| n , A 2 , m 2 -J 



(19) 



and that due to gluino exchange (Fig. lh) is 

2a 



3 7T 



TD I 2 2 2\ , / 2 2 2\d^2 2 2\ 

B (m qn ,m g , m q ) + (m ? - n - m, - m § ) B (m Sn ,m~ g , m q ) 
— 2m q mg(—l) n sm29a B n (mi , m~, m 2 ^ 



7 q I >0\," t q n i" t g>" l qJ 



(5)^2 n 



^(m^) = i—m g m, cos 20$ .B (m? n , m?,m 2 ). 



The four-squark interaction (Fig. li) gives 

4l ) «) = 4! ) «) = ^sin4^ 



(20) 
(21) 

(22) 



where A (m 2 ) is the standard one-point function in the convention of Note that S^£,(m 



sS(m|). 



91' 



3.3 On-shell renormalization of squark mixing angles 

It is necessary to renormalize the squark mixing angles 9 q by adding appropriate counterterms 
to obtain ultraviolet finite decay widths: 



Sc. 



(c) 
21Z 



-^II L cos29 q 69 q , 



(23) 



5c 



(c) 
ijW 



1 

V2 



- sin 9 q cos 9 q i sin 9 q sin 9 q / 
cos 9 q cos Op cos 9 q sin 6^/ 



+ 



cos 9q sin 6^/ 
sin 9 q sin 0$/ 



cos 9 q cos 
sin 0g cos 9qi 



59 q , 



(24) 
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For the definition of the on-shell mixing angles 9q we follow the procedure given in Ref. 0. The 
counterterm 59q is fixed such that it cancels the off-diagonal part of the squark wave function 



corrections to the cross section of e + e — > q x q 2 - 59 q 

sin 49 a 



56 f + 



59f = — 



Git m\ - ra| 2 



, ( is then given by 



59 



(.</) 



m-grriq 



TD ! 2 2 2\ ~ 

BQ{m h ,m- g ,m q )v n 



D I 2 2 2\ ~ 

B (m q - 1 ,m s ,m q )v 2 2 



(25) 
(26) 



with = 4(/3 L cos 2 ^q — e g sin 2 w ) and v 22 = 4(I$ L sm 2 9g — e^sin 2 ^). Obviously, for the 
decay q 2 — > QiZ°, the counterterm (|23D completely cancels the off-diagonal wave function 
corrections (Figs, lh (i ^ n) and li (i ^ k)). In case of q t — > q'j W the contribution of the 
squark loop, Fig. li, is cancelled. Thus, in both cases the total squark loop contribution to the 
correction is zero, 5Y® = 0. 



3.4 Real gluon emission 



In order to cancel the infrared divergence we include real (hard and soft) gluon emission: 



5T real = r(gf ^qjVg) (Figs, lj, lk, and It). The width is 



5T 



real 



9 2 c|y a s 
37r 2 m,- 



21- 



K 



I + h + m 2 fa + rrtj hi + (m 2 + m 2 - m^) J, 



01 



(27) 



Again, k = n(mf, m 2 , m 2 /). The phase space integrals J, J„, and J nm have (m«, m^, my) as 
arguments. Their explicit forms are given in [1^ . 



We have checked explicitely that the corrected decay width, eq. (]TTJ), is ultraviolet and infrared 
finite. 



4 Numerical results and discussion 



In general, the stop and sbottom sectors are determined by the soft SUSY breaking parameters 
(Mq, Mfj, and M^), the trilinear couplings (A t and A b ), fi, and tan/3, which all enter the squark 
mass matrices. In order to show the importance of the q 2 — > q x Z Q and q { — > q'jW ± decays, we 
plot in Fig. 2 the branching ratios of these modes as a function of fi for Mq = 500 GeV, 
Mjj = 444 GeV, M b = 556 GeV, and A t = A b = 500 GeV. For the SU(2) gaugino mass we take 
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M = 200 GeV, and for the mass of the pseudoscalar Higgs tua = 200 GeV (for the U(l) gaugino 
mass M' and the gluino mass we use the GUT relations M' = | Mtan 2 9 W and m g = a s /a 2 M). 
Fig. 2a (2b) is for tan/5 = 2 (40). We see that the squark decays into vector bosons can have 
very large branching ratios under the conditions (i), (ii), and (iii) given in the introduction. 

We now turn to the numerical analysis of the 0(a s ) SUSY-QCD corrected decay widths. As 
the squark couplings to vector bosons depend only on the squark mixing angles, we use the 
on-shell squark masses 2 and mixing angles 9q (0 < 9q < ir) as input parameters. Moreover, 
we take m t = 175 GeV, m b = 5 GeV, m z = 91.2 GeV, sin 2 6 W = 0.23, a(m z ) = 1/128.87, and 
ot s (m z ) = 0.12. For the running of a s we use a s (Q 2 ) = 127r/[(33 — 2n/) ln(Q 2 /A 2 t/ )] with rif 
the number of flavors. 

We first discuss the decay t 2 — > tiZ°. Figure 3 shows the tree-level and the 0(ct s ) SUSY-QCD 
corrected widths of this decay as a function of the lighter stop mass , for m 4 ~ 2 = 650 GeV, 
cos6t = —0.6, and m~ g = 500 GeV. SUSY-QCD corrections reduce the tree-level width by 
— 11.7% to —6.8% in the range of = 80 to 508 GeV. It is interesting to note that the 
gluonic correction decreases quickly with increasing while the correction due to gluino ex- 
change varies only little with m h . In our example, 5Y^> /T° = -4.5% and 6F®/T° = -7.2% 
at m h = 80 GeV, whereas at m h = 550 GeV ST® /T° ~ 0% and 8T® /T = -6.8%. 
The dependence on the stop mixing angle is shown in Fig. 4. Fig. 4a shows the tree-level width 
together with the 0(a s ) corrected width of t 2 — ► t±Z° as a function of cos# t ~ for = 200 GeV, 
m h = 650 GeV, and m g = 500 GeV. Assuming My < Mq, as suggested by SUSY-GUT, the 
stop mixing angle is varied in the range — ^ < cos# t ~ < With the tyty-Z coupling propor- 
tional to sin 26 1 the decay width has maxima at cos^ = ±-^j (maximal mixing) and vanishes 
at cosflf = 0. In Fig. 4b we plot the relative correction 5T/T for m i{i 2} = {200,650} GeV 
and m g = 500 and 1000 GeV. As the gluonic correction has the same 6~ t dependence as the 
tree-level width, SV^/T = —3% in our example, the 6~ t dependence in Fig. 4b comes only from 
the correction due to gluino exchange. For m g = 500 GeV and cos% < — 0.1 (cosdf > 0.1) the 
correction is —10% to —12% (—3% to —5%); for m g — 1 TeV and | cos^ t ~| > 0.1 the correction 
is —2% to —5%. Approaching cos^- = 0, where t 1( ^ = t R ^, 5T /T° diverges due to the vanish- 
ing tree-level coupling c 2 iz while 5c 2 i Z ^ 0. In this case the decay width becomes of 0(a 2 s ). 
Note, however, that the appearance of this divergence, as well as the condition cos# t ~ = 0, is 
renormalization scheme dependent. 

Taking a closer look on the gluino mass dependence we find that the gluino decouples slowly; e.g. 
for = 200 GeV, m t ~ 2 = 650 GeV, and cos# t ~ = —0.6 the gluino contribution to the correction 
§Y('9)/T° = -2.6%, -0.8%, -0.4% for m g = 600 GeV, 1 TeV, 1.5 TeV, respectively. On the 
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other hand, the size of the gluino contribution quickly increases when approaching the t 2 — > tg 
threshold: For m g = 480 GeV 5T® /T° = -17.2%, and for m g = 476 GeV 5T®/T° = -38.8%. 
For m-g < mf 2 — m t the gluino correction becomes positive. In our example it reaches the 
maximum at m g = 270 GeV with 5T^ /T° = 3.4%. In general, the dependence on the gluino 
mass near the threshold is less pronounced for cos# t ~ > 0. 

The decay b 2 — > biZ° can be important for large tan/3 (see Fig. 2b). The SUSY-QCD correc- 
tions to this decay behave similarly to those to t% — > t\Z°. However, the corrections due to 
gluino exchange are smaller because of the smaller bottom quark mass and thus the dependece 
of <5r/r° on the sbottom mixing angle is very weak. 

Let us now turn to the squark decays into W bosons. Here we discuss two special cases: 

(i) b\ and 62 decaying into a relatively light t\ plus W~ for small sbottom mixing (small tan f3 
scenario). In this case the mass difference of b\ and b 2 is expected to be rather small and thus 
the decays b 2 — > bi (Z°, h°, H°, A ) should be kinematically suppressed or even forbidden. 

(ii) A heavy i 2 decaying into a relatively light b\ plus W + for large sbottom mixing (large tan (3 
scenario) . 

Note, however, that in a combined treatment of both the stop and the sbottom sectors there is 
a constraint among the parameters m^, m 4 ~ 2 , % mi , m~ b2 , and 6~ b (for a given value of tan (3) due 
to SU(2) L gauge symmetry [10, 18]. Therefore one of these parameters is fixed by the others. 
Here one has to take into account that in the on-shell scheme M^(t) 7^ M%(b) at 0(a s ), see 



Rcf. 

The decay widths of 61.2 — > t{W~ are shown in Fig. 5a as a function of the stop mixing angle, for 
m ~ bi = 500 GeV, m~ h = 520 GeV, cos% = -0.9, m h = 200 GeV, m g = 520 GeV, and tan/5 = 2. 
The value of m~ t2 is determined by the other parameters as discussed above. Hence m~ t2 varies 
from 533 GeV to 733 GeV depending on the stop mixing angle in Fig. 5a. (Here we have also 
checked that the squark parameters do not cause too large A t or Ab to avoid color breaking 
minima.) Despite the larger phase space for the b 2 decay, the width of 62 - * tiW~ is smaller 
than that of bi — > t\W~ because the W couples only to the "left" components of the squarks 
(note that b± ~ b L and b 2 ~ b R for cos#^ = —0.9). Fig. 5b shows the relative correction ST/F 
for the same squark parameters as in Fig. 5a. For m g = 520 GeV and | cos^-| > 0.1 SUSY-QCD 
corrections are —11% to +4%. For m g = ITeV the effects are weaker, i.e. about —5% to — 1%. 
Again, there is a "singularity" at cos^ = where the tree-level t 1 -b i -W ± coupling vanishes. 
The dependence of <5r/r° on the sbottom mixing angle is, in general, much weaker than that 
on the stop mixing angle (apart from a singularity for the b { being a pure b R ). The overall 
dependence on the gluino mass is in general similar to that of the t 2 — > t x Z® decay. However, 
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the enhancement effect of the threshold at m g = mg. — mb is less pronounced. 
An example for large tan (3 is shown in Fig. 6. In Fig. 6a we plot the tree-level and the SUSY- 
QCD corrected widths of t 2 — > b\W + as a function of cos % for = 300 GeV, m t ~ 2 = 650 GeV, 
= 380 GeV, cos^ = —0.8, m g = 500 GeV, and tan/? = 40. m~ b2 is calculated from the 
other parameters and thus varies from 615 GeV to 918 GeV. As expected, the decay width is 
maximal for t 2 = t L and vanishes for t 2 = t R . In the example chosen, SUSY-QCD corrections 
are -2.4% to -4.7% for m g = 500 GeV and about -1% to -1.5% for m g = 1 TeV as can be 
seen in Fig. 6b, where 5T/T is shown as a function of cos6> t ~ for the same squark parameters 
as given above, and m g = 500 GeV and 1 TeV. Again, there is almost no dependence on cos 9~ b 
apart from the singularity at cos9~ b = 0, i.e. b ± = b R . 

5 Summary 

We have calculated the 0(a s ) supersymmetric QCD corrections to squark decays into vector 
bosons in the on-shell renormalization scheme using dimensional reduction. In particular, we 
have discussed examples for the decays t 2 — > hZ , t 2 — > b{W + , and b lj2 — > t{W~ . We have 
found that the correction 8T/T is typically of the order —5% to —10%, depending on the squark 
masses, squark mixing angles, and m g . Near the q — > qg threshold the correction can also exceed 
—10%. It has also turned out that the gluino decouples slowly. Moreover, the gluino-exchange 
corrections alter the 9 g dependence of the tree-level widths. For squark mixing angles where 
the decay width vanishes at tree-level, the gluino corrections do not vanish and may lead to 
non-zero widths of 0(a 2 s ). 
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Fig. 1: Feynman diagrams relevant for the 0(a s ) SUSY-QCD correc- 
tions to squark decays into vector bosons. 
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Ii [GeV] 

Fig. 2: Tree-level branching ratios for stop and sbottom decays into 
vector bosons as a function of fx, for Mq = 500 GeV, Mjj = AAA GeV, 
Mfy = 556 GeV, A t = A b = 500 GeV, M = 200 GeV, and m A = 200 
GeV. (a) tan/3 = 2, and (b) tan/3 = 40. 
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100 200 300 400 500 



m lx [GeV] 

Fig. 3: Tree-level (dashed line) and SUSY-QCD corrected (solid line) 
widths of the decay t 2 — > tiZ° as a function of for m^ 2 = 650 GeV, 
cos &i = —0.6, and rrig = 500 GeV. 
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Fig. 4a: Tree-level (dashed line) and SUSY-QCD corrected (solid line) 
widths of the decay t 2 — > tiZ° as a function of cos Q~ t for = 200 GeV, 
m h = 650 GeV, and m g = 500 GeV. 
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Fig. 4b: SUSY-QCD correction relative to the tree-level width of the 
decay t 2 — > i\Z° as a function of cos^- for = 200 GeV and m t - 2 = 650 
GeV. The black (gray) lines are for m- g = 500 (1000) GeV. 
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Fig. 5a: Tree-level (dashed^ lines) and SUSY-QCD corrected (solid 
lines) widths of the decays &i )2 — > tiW~ as a function of cos% for 
m 5i = 500 GeV, m 52 = 520 GeV, cosflj = -0.9, m fl = 200 GeV, 
m§ = 520 GeV, and tan/5 = 2. mf 2 is a function of the other parameters 
and varies with cos#f. 




Fig. 5b: SUSY-QCD corrections relative to the tree-level widths of the 
decays bi — > t{W~ (dashdotted lines) and 62 — ► t{W~ (solid lines) as a 
function of cos#£ for = 500 GeV, = 520 GeV, cos^ = —0.9, 
ra^ = 200 GeV, and tan/3 = 2. The black (gray) lines are for m~ g = 520 
(1000) GeV. m t ~ 2 is a function of the other parameters, as in Fig. 5a. 
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Fig. 6a: Tree-level (dashed line) and SUSY-QCD corrected (solid line) 
widths of the decay t 2 — > biW + as a function of cos Of for = 300 GeV, 
m h = 650 GeV, m bi = 380 GeV, cos0 b = -0.8, m g = 500 GeV, and 
tan j3 = 40. m b2 is a function of the other parameters and varies with 
cos 6~ t . 
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Fig. 6b: SUSY-QCD corrections relative to the tree-level width of 
the decay t 2 — > biW + as a function of cos% for = 300 GeV, 
m h = 650 GeV, m~ bi = 380 GeV, cosflg = -0.8, and tan/3 = 40. The 
black (gray) line is for m~ g = 500 (1000) GeV. m~ b2 is a function of the 
other parameters, as in Fig. 6a. 
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